Linear theory of unstable growth on rough surfaces 
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Unstable homoepitaxy on rough substrates is treated 
within a linear continuum theory. The time dependence of the 
surface width W(t) is governed by three length scales: The 
characteristic scale lo of the substrate roughness, the terrace 
size Id and the Ehrlich-Schwoebel length Ies. If Ies *C Id 
(weak step edge barriers) and lo <C lm ~ Id\/Id/Ies, then 
W(t) displays a minimum at a coverage # min ~ (Id/Ies) 2 , 
where the initial surface width is reduced by a factor lo/l m - 
The role of deposition and diffusion noise is analyzed. The 
results are applied to recent experiments on the growth of 
InAs buffer layers [M.F. Gyure et al., Phys. Rev. Lett. 81, 
4931 (1998)]. The overall features of the observed roughness 
evolution are captured by the linear theory, but the detailed 
time dependence shows distinct deviations which suggest a 
significant influence of nonlinearities. 



I. INTRODUCTION 

A high symmetry crystal surface growing epitaxially 
from a molecular beam can become unstable towards 
the formation of pyramidal mounds if the mass trans- 
port between different atomic layers ia, reduced by ad- 
ditional energy barriers at step edgesoa. Over the last 
few years, this phenomenon has been observed for a wide 
range of metal and semiconductor surfaces, and a consid- 
erable body of theoretical work has been devoted to the 
description of the asymptotic (late time) evolution of the 
surface morphologyuEl. In the early time regime contin- 
uum theory predicts an exponential growth of the surface 
modulations. For this reason the precise initial state of 
the surface has commonly been disregarded, since the ex- 
ponential instability should rapidly wash out the details 
of the substrate roughness. 

In a recent papeiQ Gyure, Zinck, Ratsch and Vveden- 
sky (GZRV) presented experimental and numerical re- 
sults for the early time development of unstable homoepi- 
taxy from a rough substrate, which show a more complex 
scenario: It was observed that the competition between 
smoothening of the initial roughness and the instability 
associated with the incipient mound structure can lead 
to a minimum in the total surface width. A similar effect 
was predicted previously in the context of noise-induced 
roughcningu, and related experimentaLpbservations have 
been reported both for thin metal filmsa and semiconduc- 
tor multilayersa. Qualitatively, the minimum originates 
from the wavelength dependence of smoothing and (de- 



terministic or stochastic) roughening rates: If the rough- 
ness spectrum of the substrate has sufficient weight at 
short wavelengths— which are efficiently smoothened by 
capillarity effects^], then the decrease of the substrate 
contribution to the surface width can temporarily domi- 
nate the long wavelength roughening induced by growth. 
The possibility to minimize the surface roughness by 
an appropriate choice of the buffer layer thickness and 
other growth parameters is of obvious interest in appli- 
cations. In this paper we develop a quantitative the- 
ory of unstable growth on rough substrates, which allows 
us to determine the conditions under which a minimum 
occurs, and to estimate the layer thickness of minimal 
roughness in terms of microscopic length scales and pa- 
rameters, such as the in-layer and interlayer diffusion bar- 
riers. Our starting point is the observation of GZRV that 
the time evolution of the roughness spectrum appears to 
be well described by the linearized continuum evolution 
equatiop_for the surface. By incorporating various kinds 
of noiseEa into the linear theory we can compare the in- 
fluence of stochastic and deterministic roughening, and 
obtain a unified description of both cases. A critical dis- 
cussion of our results in relation to the experiments of 
GZRV will be presented at the end of the paper. 



II. LINEARIZED CONTINUUM THEORY 



The standard phenomenological evolution equati 
the continuous surface profile H(r,t) is of the form! 



d t H + V • J = F 



(1) 



s both a growth- 
2 ' and a stabilizing 



where the surface current J incor 
induced destabilizing contributio: 

term originating in capillaritylMa, and F denotes the depo- 
sition flux, which will be assumed constant for the time 
being. Small fluctuations h(r, t) around the flat singular 
surface H — Ft then satisfy the linear equation 



d t h = -aV 2 h - K{V 2 fh 



(2) 



with positive coefficients a, k representing deposition (a) 
and smoothening (k), respectively, whose relation to the 
growth parameters will be explained below. 

The substrate roughness is incorporated through a spa- 
tial roughness spectrum (\h(k, 0)| 2 ) = S(k,0) = So(k), 
where h(k, t) is the Fourier transform of h(r, t) and 
k — |k|. Under the linear equation (||) the roughness 
spectrum evolves as 



S(k,t) = S {k) exp[2(afc 2 - K k 4 )t] 



(3) 



which implies that fluctuations with wavenumbers k > 
k c = yj aj k arc damped while those with k < k c are 
exponentially amplified. The surface width W(t) is ob- 
tained by summing over all wavenumbers, 



W 2 {t) = 2ir dk kS (k)e 
Jo 



2(ak 2 -Kk i )t 



(4) 



Motivated by the experimental data shown in Figure 4 
of GZRV we choose a white noise roughness spectrum, 



*w = { «"*/■* 



k < ir/l 
else, 



(5) 



where the small scale cutoff Iq is required for a finite value 
Wo = W(0) of the initial surface width. Taking the time 
derivative of (Q) and evaluating it at t = 0, we find that 
the surface width shows an initial decrease if U, m /lo) 2 > 



12, where l m — 2-k^/2kJol is the wavelengthQ of those 
fluctuations which are maximally amplified by the linear 
equation (0). Thus the condition for a nonmonotonic 
time dependence of the surface width is that the length 
scale characterizing the substrate roughness, Iq, is much 
smaller than the typical scale l m of the emerging mounds. 
This result holds also for more general initial roughness 
spectra, e.g. So(fc) = Ak~ p with p < 2 and a small scale 
cutoff Iq. For substrates whose roughness is dominated by 
long wavelength fluctuations, in the sense that So ~ k~ p 
with p > 2, a large scale cutoff is needed and the time 
derivative dW/dt\ t=0 does usually not exist. 

In the following we take Iq <§C l m - Then Eq.(0) reduces 
to the scaling form 



W 2 (t) = W 2 (l /I,n) 2 ^(t/1 



(6) 



where 1/r = a 2 /4k is the amplification rate of the max- 
imally unstable fluctuations and the scaling function is 



$(*) = e 2x y/2w/x[l + erf(v^)] 



(7) 



with erf (s) = (2/^/tt) L exp(— t 2 )dt. The width attains 
its minimum at a time i m i n « 0.18 r, where it has been 
reduced by a factor 



W(t min )/W Q « 3.7 (lo/lrn). 



(8) 



Since the factor 1 + erf (y/2x) in (0) only varies between 
1 and 2, the scaling function $(#) is essentially the prod- 
uct of a decaying power and an exponentially increasing 
factor. The power law for small x reflects the particular 
smoothening mechanism (capillarity-driven surface diffu- 
sion) and its general formEl is given by ([l8|) below. For 
finite lo/lm the power law sets in for times t > to with 
to « (lo/lm^r. 

To relate the behavior of W(t) to microscopic parame- 
ters we need to express the coefficients a and k of (g) 
in terms of the two length scales governing unstable 



homoepitaxyerEl: The typical terrace sizet2l Id and the 
Ehrlich-Schwoebel-lengthE3 



l ES = a,, (£>/£>' - 1) = a||(e A£ / fc - T - 1) 



(9) 



defined in terms of the in-layer lattice spacing an, the 
in-layer (interlayer) surface diffusion constant D (D 1 ) 
and the step edge barrier AE. Comparison of the two 
length scales allows to distinguish conditions of strong 
{Ies ^> Id) and weak (Ies *C Id) step edge barriers; 
in the first case a « Fl 2 D , in the second a rj FIdIes- 
The coefficient k is traditicuially associated with near- 
equilibrium surface diffusionlij, however under far-from- 
equilibrium growth conditions the dominant contribution 
to k is— believed to arise from the random nucleation 
processtH The expression k ~ Flf-> is then suggested 
by dimensional analysisEJ aad scaling argumentalj. It 
leads to a consistent pictures in the sense that l m w Id 
and r w F^ 1 in the strong barrier case, which implies 
that mounds develop on the submonolayer islands al- 
ready during, the .growth of the first few layers ( "wedding 
cake" regimeEjii3) . In the weak barrier case we find 



/, 



IdVId/Ies and T~F- x {l D /l ES ? 



(10) 



The minimum in the surface width thus occurs at a cov- 
erage 



(Id/I 



ES) 



>1 



(11) 



which corresponds, not surprisingly, to the coverage 
where mounds ficst, become visible for growth from a 
smooth substrateatJ. Similarly the coverage 9q = Fto 
at which the scaling form (|^) for the width begins to 
hold is of the order of 8 ~ {Io/Id) 4 independent of Ies 
(provided l ES <C Id)- 

To apply these considerations to the experiment on 
InAs growth of GZRV, we first need to check the condi- 
tion Iq <§C l m . ^From Figure 4 of the paperQ we estimate 
that VFo/W^limin) ~ 4. Comparing this to the theoretical 
prediction (||) we find Iq rs 0.07 x l m and l m /lo ^ 1 is 
true. This is in contrast to the kinetic Monte Carlo simu- 
lations of GZRV, where Wo/W(t min ) w 1.1. The instabil- 
ity length in the experiment is l m ~ 1.0 pm, which yields 
lo ~ 70 nm for the small scale cutoff of the substrate 
roughness. This is consistent with the initial roughness 
spectrum in Figure 4 of GZRV, which is constant at least 
down to a length scale of 300 nm. The minimum width 
is attained at a film thickness of about 0.57 pm. Using 
aii ~ 6A and a bilayer thickness a± fa 3A, we therefore 
estimate that 9 m [ n « 1900 and l m /a\\ ~ 1700, and hence 
Ies/c^w ~ 6 and Id/o,u ~ 250. At the experimental tem- 
perature of 500° C this implies a step edge barriex_Ai? 
of the order of 0.1 eV, comparable to estimatesSO for 
GaAs. 



III. NOISE EFFECTS 

Next we include a noise term r](r,t) in Eq. (ph. The 
different sources of noise, the individual events of depo- 
sition ( "shot noise" ) and diffusipfo-enter the noise corre- 
lator with different dependenceOiIa on the wavenumber 
k. We write it in the form 



R(k) = (tj(k,t)»?(-M)) = Rs + Rvk 2 



(12) 



with Rs and Rd denoting the strength of deposition and 
diffusion noise, respectively. In the linear model with 
noise the roughness spectrum S(k, t) then contains a part 
reflecting the history of the noise, as well as the deter- 
ministic evolution of the initial roughness treated above. 
The full expression reads 



S(k,t) = Sdct(M) +Snoisc(M) = 



(13) 



5*o (k) e 



2{ak 2 -Kk i )t 



R(k) 



2(ak 2 - Kfc 4 ) 



D 2(afc 2 -refc 4 )i 



Unlike the deterministic mechanism, the noise increases 
the amplitude of the spectrum for every wavelength, i.e. 
dtS no [ so (k, t) > for all k. We shall now examine whether 
under the experimental conditions of GZRV noise sub- 
stantially contributes to the surface width. 

The deposited particle flux can be seen as a Poisson 
process with intensity F, so Rg — a±a 2 F. Shot noise 
thus contributes to the total width by 



W s (t) 2 = FTa x {au/l m ) 2 *(tM, 



(14) 



with 4»'(x) = $(») for the choice (§) of S (k). Using M) 
we see that (|l4]) can be can be ignored against Eq. M) if 



(Wo/a ± ) 2 (Z /a||) 2 » (Id/I 



ES) 



(15) 



With our estimate Iq « 70 nm this condition is satisfied 
in the experiment. A different interpretation of Eq.fllq) 
will be given below. 

The diffusion noise strength is guzen by the average 
rate of adatom jumps on the surfaceO, so 



R D ~ Pl D~l 2 D F 



(16) 



where we have |Used the estimate p\ ~ Fl D /D for the 
adatom densityu. Diffusion noise thus becomes more 
important than shot noise for k > it /Id, whereas it can 
be neglected for long wavelengths. For large k we can 
approximate the contribution of diffusion noise in ( |13] ) 
by Ri)k 2 / (nk 4 ) which enters the total width as Wr>(t) - 
l D /ifnFr \og(lo / a\\) . This means that roughly Wd(£) 2 ~ 
(lD/l m ) 2 Ws(T) 2 <§; Ws(t) 2 because l m > Id in the weak 
barrier regime. In particular, at the time when the width 
minimum is attained, diffusion noise can be neglected 
against shot noise and for the experiment of GZRV Eq. 
(a) remains valid. 

It was mentioned already that due to noise a minimum 
in the surface width may occur even in the absence of 



step edge barriersQ. For completeness we provide here 
a simple analysis for the most general linear Langevin 
equation of kinetic roughening, 



d t h 



-{-vV 2 ) z/2 h + rj 



(17) 



where z — 2 and z — 4 correspond to evaporation- 
condensatioJ3-|and surface diffusion dominated relaxation, 
respectivelyli3, v > is a constant and and r\ is the de- 
position noise. Odd or noninteger values of z describe 
nonlocal relaxation mechanisms and can be treated on 
the same footingQ. The linearity of ( p7| ) implies that the 
substrate contribution and the growth induced contribu- 
tion to the roughness can be separatedD. The substrate 
contribution is found to decay according to 



WW*) « w (io/m) d/2 



(18) 



for a d-dimensional surface, for times such that the cor- 
relation length of the growth-induced roughness £(£) ex- 
ceeds Iq (otherwise W su b ~ Wo)- In terms of physk 
cal quantities the correlation length can be written astil 
£,(t) ~ IdQ 1 ^ ■ To determine the coverage m ; n of mini- 
mal surface width for purely stochastic roughening, the 
substrate contribution (fy|X should be compared to the 
growth induced roughnesscJ 



TU 



rowth 



a ±(e/e)( z - d)/2z 



(19) 



where 8 is the coverage at which the width becomes of 
the order of a± and thus lattice effects (such as tem- 
poral oscillations of the step density) die out; the_ex^ 
pression (M) holds for 6 3> 9. With the estimatai3'fL2l 
6~ {l D /a\\) zd ^ z -^ we obtain 



e„ 



(W /a ± ) 2 (l a / ail ) d 



(20) 



independent of z. Comparing Eqs. (EOT) and (O) we 
have thus recovered the crossover condition ( |l5| ) be- 
tween deterministic instability and stochastic roughening 
from the opposite side. To neglect the growth instabil- 
ity in Eq.(O) is no longer justified when the minimum 
width coverage 9 m i n predicted by the deterministic the- 
ory (Eq.(0)) is smaller than 9 min in Eq.(po|). 



IV. DISCUSSION AND CONCLUSION 

The main results of this paper are Eqs. (g) and J14), 
which express the time dependent surface roughness in 
terms of the characteristic length and time scales of the 
problem - the substrate roughness scale Iq, the incipient 
mound size l m , and the linear growth time r -, the lat- 
ter two of which are, in turn, related to the microscopic 
growth parameters through Eq. (UG) . For the experiment 
of GZRV, the measured values of l m and r were seen to 
imply reasonable numbers for the microscopic lengths Id 
and Ies, an d for the step edge barrier AE. 



It is then natural to ask to what extent the linear the- 
ory can be used to quantitatively describe the experi- 
mentally observed roughness evolution, beyond provid- 
ing consistent order-of-magnitude estimates. Inspection 
of the data for the surface roughness depicted in the inset 
of Figure 4 of GZRV quickly leads to the conclusion that, 
despite a similar overall appearance, the shape of W(t) is 
not well reproduced by our scaling functions (|6|) and ( 14 ) . 
In fact the data for the structure factor in Figure 4 show 
a qualitative feature which the linear theory is unable to 
explain: It is an immediate consequence of Eq. (O) that 
S(k, t) is a monotonic function oft (increasing or decreas- 
ing) for any k; in contrast, the measured structure factor 
shows a nonmonotonic dependence on film thickness for 

This prompts the question whether the use of the lin- 
earized theory is really justified under the experimental 
conditions. Nonlinear-terms in the surface current J in 
Eq. ([j]) are expectedQ'ErE] to matter when the surface slope 
\Vh\ becomes comparable to o,±/Id- Since typical slope 
values of the initial surface profile are of the order of 
W /l , the condition for the validity of the linear theory 



W /aj_ < Io/Id- 



(21) 



In the experiment of GZRV Wo/a± ~ 3 and, from the es- 
timates presented above, Io/Id ~ 0.5; thus the condition 
( p2l| ) is Uveakly) violated. The analogy to phase ordering 
kincticsu suggests that the early time evolution may Jae 
qualitatively altered when nonlinearities are importantE! 
This seems like an interesting problem for future study. 
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